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ABST RAG I 


A new method for predicting the off design performance ofa 
multistage compressor has been developed. This method replaces 
the usual routine numerical integration through the compressor by 
approximate analytical integration. The overall compressor perform- 
ance of the compressor is found as a perturbation of the design 
performance in terms of linear and quadratic expansions in two 
parameters, the deviation from the design flow rate and the deviation 
from the design speed. The procedure for this preliminary examina- 
tion has been simplified by assuming that all the single stage 
performance characteristics are identical and that the difference 
between static and stagnation conditions are negligible. 

The results of this method are quadratic equations,in terms of 
the two perturbation parameters,for the temperature ratio, pressure 
ratio and adiabatic efficiency of the overall compressor. The 
amount of numerical work in finding the coefficients of these quadratic 
expansions is not negligible, however the method appears promising. 
Reasonably accurate results are expected for speeds varying as much 
as ten per cent from design and flow rates varying as much as four 


per cent from the design flow. 
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I, INTRODUCTION 


The usual technique for calculation of multistage compressor 
performance is numerical. Starting at the entrance with a prescribed 
flow rate and rotor speed the pressures, temperatures and flow 
coefficients are determined behind each stage in succession from the 
known performance characteristics of the individual stages. Although 
this method can be carried out rather rapidly on machines, it is 
desirable to have available a method that is more suitable for explor- 
atory investigations of the influence of stage characteristics. 

A new method is proposed which replaces the routine numerical 
integration through the compressor by approximate analytical integra- 
tion, The overall performance of the compressor can be found asa 
perturbation of the design performance in terms of linear and quadratic 
BE ansions in two parameters, the deviation from the design flow rate 
and the deviation from the design speed. The amount of numerical 
work in finding the coefficients in the expansion is not negligible, 
nevertheless the method appears promising. 

The investigation should be considered as preliminary. Cer- 
tain approximations, such as the assumptions that all single stage 
performance characteristics are identical and that the difference be- 
tween stagnation and static conditions are negligible, are not essential 


but simplify the procedure for the initial examination. 
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HW. DEVELOPMENT OF EQUATIONS 


~~ Work Coefficient 

The performance of multistage compressors is generally 
presented as plots of temperature and pressure ratios against dimen- 
sionless weight flow for given dimensionless compressor speed. The 
stage performance is most conveniently represented in the dimension- 
less form as a work coefficient ('), pressure coefficient (W’) and 
adiabatic efficiency (yn) against a flow coefficient (o). These 


dimensionless quantities are defined as: 


vy - Ah. _ cpAT 
U5 


y’ = 
p 





 |[> 
ab SS 


mM 


OUA 


my = 


<8 |x 


Ah. is the increase in stagnation enthalpy per stage, s the specific 
heat of the gas at constant pressure, AT the increase in stagnation 
temperature per stage, Ap the increase in stagnation pressure per 
stage, and p the gas density. The compressor rotor reference 
speed is represented by U and the mean axial gas velocity by A i, 
m is the mass flow of gas and A the cross sectional area of the 


flow annulus. 








3 
The multistage compressor temperature ratio, pressure ratio, 
and adiabatic efficiency will be determined from the single stage 
dimensionless characteristics. 
In general it can be said that the work coefficient Y isa 


function of the flow coefficient: 


Y = FW) 


The temperature rise (AT) ina single stage can thus be writtenasa 


function of the flow coefficient: 


Eo UE URC 
Oi tae Te F (9) 








A multistage compressor will be assumed to be made up of 
stages which all have the same characteristic. Such a compressor 
can be considered to be made up of an infinite number of infinitesimal 
stages such that the temperature can be considered to be a continuous 
variable along the axis of the compressor. Let the continuous axial 
variable be €, replacing the sequence En = y where n is the 
stage number from the entrance and N is the total number of stages 
of the compressor (0 £n<N). Using the variable € the tempera- 


ture rise can be written in the differential form: 


dT u* 
a = YO rg) 
dé Cp 
The compressor then has an overall "length" of unity, and the 


function ‘TT (6) becomes: 


= F (9) NW 
T(>) = N . : 





4 
The adiabatic compressor efficiency is normally defined as: 


ia WES 
ii Re 


where T is the actual temperature and T 


2 


i is the isentropic final 


temperature ie a given avin rise. Introducing the isentropic 


relationship ae (eye this becomes: 
, (Ys) 7 

T/ Te) — | 

here © and T/T(o) are the total pressure and temperature 
Ww e ae (AG) eer pr ie mp r 
ratios, p(o) and T(o) are the initial pressure and temperature and 

¥ is the ratio of specific heats of the gas. Ina single stage the 
pressure and temperature increase can be represented by Ap and 
AT suchthat p = p(o)+ Ap and T = T(o)+ AT. The stage 


efficiency thus becomes: 


(| i AP/o_) F — 


VAS AG 
T (9) 





% = 





Since Ap is small compared to p(o) the pressure term can be 


expanded to give: 





Therefore in the limiting case, with Ap very small (as in the single 


stage): 





yo i AP/p 0 


A T/T) 


or 
i ae i Al 
Ae y-|/ TC) a 
This will also be the form of the limiting case as Ap-—»o so that the 


differential form can be written as: 
i dp =— 7 my. a _dT (Tay 
ea & a) © le 

The adiabatic efficiency (yn) as defined here can be shown to 


be equivalent to the definition » = /ay in the limiting case. By 


the definitions of the work and pressure coefficients: 


ae 
wo ey) AP 
VY Coeee ae. Aan 
ie CSP 
But for Ap small equation (1) applies and: 
y )& 
i iia ( y-l/ TT 
¥ C oP 
Substituting the thermodynamic relations for a perfect gas: 
eT Cp y 
gives j 
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B, Differential Equations for Flow Coefficient 
The flow coefficient has been defined as: 
b = aol 
- e UA 
Introducing the perfect gas relation, 


p=eRt 


the flow coefficient becomes: 
-, al 
p pUuA 
Taking the logarithm and differentiating with respect to € this 


becomes: 


ep 
| 
iz 
ae 


. Eo ee 


where the ( ~) denotes differentiation with respect to the variable €. 


: 7 
Substituting for " from equation (la) this reduces to: 


eee eee ee. ES 
: jE ; (2) 


This is an approximate relation since the differences between static 





and total temperature and pressure have been neglected in the 
development. 

At the compressor design point the flow coefficient, efficiency 
and temperature rise are assumed to be constants throughout €, and 
equal to the single stage parameters at the design point. Let the sub- 
script ( ,) denote the design condition. Therefore at design 


conditions: 








bd = %, = constant 


n = No = constant 


(aT) = T= Us + constant 


(3) 





At the design point equation (3) can easily be integrated to give the 


temperature at any point in the compressor as: 


rs 
ee = (4) 


oP 
where ®, is the inlet temperature. 


Dividing equation (3) by (4) gives: 





| Uy a 
: i a = 
—- § - © (> 
where 8 is defined as: 
ee = nel Se 
° Ut. 


Equation (2) can now be solved to determine the area variation 


of the compressor: 


A = — {26% -i} a _o, 
A fo Te Q. 
Letting the term no) -~l = f£ and noting that b. = 0 (since >. 


is a constant) and substituting in equation (5): 


A ie (6) 


Smee ee 


A caer: 





8 
This area distribution is fixed by the design geometry. Substituting 
(6) into equation (2) the general differential equation becomes: 
nb Opin. pi Fr (7) 
> T E+ 6. 


To eliminate the temperature (T) equation (7) is solved for T and 


substituted in the expression for T: 














am 
l= f, 
Ete, e 
; c. 1 
fe cs ea OO, See zi _ dj_t ge 
U jof= dé fo eos die f. _ oF 
E+@ oo E+B. 
~- ae lee 
gE ome wae kOe 
E+ 6. b 


Performing the differentiations and substituting the relations: 





ies dens ow 
dE dd dé do 
he, he el le by 
: do dé a6 


there results the general differential equation in 6: 
Ep T) f, ae OF a a7) 2] 7 
uot ages (rs larry au ‘ 


ler + fap) e cis + rd = O 


(8) 
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C. Solution for the Flow Coefficient 
A series expansion of each of the terms involved in the 
differential equation has been used to find a solution for the flow 
coefficient. From the single stage performance it is seen that both 
the efficiency function f and work coefficient 7 are functions of 
the variable 9: 
f= F(¢) + = T($) 
These functions can thus be expressed as a Taylor series expansion 
about the design point. The stage efficiency at designis assumed to 
have its maximum value, hence the first derivative of the efficiency 
function becomes zero ey) = 0. The Taylor series expansions 


do 


thus become: 








{(>) = f |. +(44) (¢ _$) 4 + (44) (o-a)+ eee 


(9) 
rele T+ (SE)CO de 2 (ED \o-ge 2 (ft )@-ahe ~~ 


The flow coefficient will be represented by a series of the 


form: 


bo = >, +(6->)— db + epred +6$ + Soe CD) 


The factor © represents a convenient small parameter to be 
replaced by unity in the final form. Substitution of the series (9) and 
(10) into equation (8) results in a series equation in powers of €. 


Equating terms in like powers of € gives a series of linear 
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differential equations in the flow coefficients >, (ni =o), saa 


= ie 
Only the first two of these coefficients will be retained. The linear 


differential equations for these coefficients are: 








>, + {f a? +2] ze, 9,7 0 (11) 


b=] 





> + 4h 4, +t b= 1e 4, i 126} @ 








+ B(e4+2 - 42 —4(B+ $14, - 4 


{ “ j 2 
7 Fs +B, P. 


12 (12) 


Note that the ({ . denotes differentiation with respect to the varible 
: when used with the flow function and denotes differentiation with 
respect to the variable ¢ when used with f or T. The constants 

5? ie ’ iF , fo fe and b, are all known from the single stage 


characteristic. 


To simplify the solution of equations (11) and (12) the following 


notations have been introduced: 





_ Tres 
52). ae 








1 


The coefficients of the left side of equation (12) become: 





B= e{at+htl] 5 CH Et) -i1¢; Dee 


Using this notation equations (11) and (12) become: 





uw ror / 
+ E+ he aes 


: oO fe b, >, (12a) 
K+ Fret Bh+ C E+ pe +D ghar 





ot 
Multiplying equation (lla) by the integrating factor (€ + B,) 


(ge +p) 0+ «(6+ A) $'— 0 


4-{ (g+a) | we 


Integrating: 
a , 
fort (eb ) > = constant 


The constant can be found by introducing the boundary condition 


at € = QO; >, = >) (o) and therefore: 


Constant = Bo b, Co) 


Therefore: 


% = (a + Bo ) $<) pce) 
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Integrating between the limits 0 to €: 


4. Rohs J ¢ -) ‘3 
b = (0. { _ (seh) "| d (0) t (0) 


a | ey 








In terms of J, this becomes: 





+I 
(D0 / 
) = (0) — ot (2) a i P (0) (13) 


E,quation 12a is solved ina similar manner. The integrating factor 
fe + ome, igs again applied to the left side of the equation and the 
solutions for b, and $, form (11b) and (13) substituted into the 


right side. On integrating once the equation is of the form: 


“24, + 3) -(%,+2) £42) 
' oa H In 8 Pe 
>, — E(é- 6.) F(g+p.) | +G(E+@) + (Exe,)” : 


The coefficients EF, F, G, and H area collection of all the constants 
in the equation not involved in the integration. On integrating the 
second term of the flow coefficient becomes: 


-2(6,¢+) -(%.+l) E+] 


E 2(%,+ 6 4 
= 9, (0) 7” E 20,1) { (-4 if, = y= Fz +I a 


-(f, +1) Xt l 
-G In a al = an ~ {Ge nga, - malt zai] 





(14) 














The coefficients of equation (14) are defined as follows: 





13 
2G, +2) P 2 


0 
= {B- eee TEN I+! ? oe 
CY Or bi? 
F = {B- se) ve ‘& a7 6 : ti 0 
Ee oe (0) + @,  b, 601 


D 2 , & : ae 
G= Garp & @ (0) + re peor > 60} + oy d ©) 











(Ts) 








Z.+3 


Ly = {c- 2 2 pl L— ber + {c - 2D} oy" * cor $ col 











All terms containing the variable & will be designated by lower case 


b's with subscripts: 








25, +1 

ie wa { (2%) : if 
ee 

oe #4 ( ra) ; | 
oF 














2(k, 
CF ° 
b= alge) 7 i} ve) 
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Introducing the coefficients (15) and (16) into equations (13) 


and (14) the first and second flow coefficient terms become: 


d (0) ee b, Q Co) (13a) 


p, (0) + q bor te po) b+ d (0) +2 dco) (14a) 


where the coefficients of (14a) are defined as: 


= £5 {lB - Gees) + bs [B+ + ane at)t byb,(C - 22) +btal 2} 











c= b(~ 5; + PP) thb(c-2 wan leas cra - 2 


dip (b, — eb) eae 


Since a large number of coefficients will be used in the 
following equations, they will be represented by lower case letters 
with number subscripts. The letter will be used to indicate coeffi- 
cients of a given term (a for $5 (0)s Cmior ; (0) >, (0); d for b* (0) 
and e for either 6; (0) or $5(0)). The number subscript refers to 
a given equation. The lower case b refers to terms containing the 


variable €. 


D. Flow Coefficient Initial Conditions 
Equations 13a and 14a can be used to determine the flow 


coefficient in terms of the initial conditions at € = 0; (, (0); $1 (0), 





Ins 
> (0) and $,(0)). The first stage of the compressor will be assumed 
to operate with the same characteristic as the single stage, therefore 
at € = 0, (0) will be known. Let (0) = $6 7 E>, (0) with (0) =0 


for n 71, thus  , (o) is prescribed and , (0) = Q. For the initial 





condition at € = 0 equation (7) becomes: 
/ To) ¢ 
== {(0) —!— (0) (17) 
Do) f To! + B, } D 


By definition: 





ea eri | eee 
qa ( dé |. me 2 ac oe 
TO= O&O 
; (18) 
Te — Ute) ut eg 
T Co) Cp cS) a CoO Te 


a 


The term ui. is similar to the term defined ash , therefore let: 
Cp & C. 





a cp © m 19 
eu: eo 4 


Substituting from the expansions (9) and (10) gives: 
dom b + € PO) {pcl—o fer n>l} 
/ 2 4 = 2 
Too) = t+ER P+ € = b (0) 
“ 2 
foo) = Ff + e* te (0) (20) 


(0) = € d “Co} + ‘ p, co) 
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Substituting (18) and the expansions (19) and (20) into equations (17) 
results in a series equation in powers of © . Equating terms of like 
powers of € gives the following expressions for the flow coefficient 


derivatives at the initial conditions: 








io) = ae 6, a d (0) 


b (ola (+i l,4 oye £. 4,00 ae 4 ra be 


lay Flow Coefficient - Final Form 

Substituting the initial conditions (21) into the expressions for 
the flowcoefficient terms (13a and 14a) results in a quadratic 
expression for the flow coefficient in terms of the velocity (or speed) 
perturbation By , and the flow coefficient perturbation >, (0) of 


the form: 


2 5 x 
(6) = >, HEALG +e ARB +e Anping EAL GOL And, 


The upper case letters will be used to represent coefficients of a 
given equation and the subscripts refer to the terms in the equation. 
The subscripts 1 and Z will be used to represent coefficients of the 
terms , (0) and 6, respectively. 

Small deviations from the design point are accomplished by 
letting € , (0) and € By be small. This can be accomplished by 


letting € = 1 and using small perturbation of >, (0) and §).-. 
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Letting € = 1 the flow coefficient becomes: 


b= bt Arh +An Or HAL GO +A, 400+ A, O00) (22) 


The coefficients of equation (22) are expressed in terms of the initial 
design conditions CRs uae fate b, ) and the location along the length 
of the compressor €. These coefficients are most easily expressed 
in terms of the coefficients of all the previous equations involved. 


Express the initial conditions (21) as: 


db wo) = cS, $ (0) + B, 6, 


(21a) 


p, (0)= 15), ,°(o) or IB G, $, (0? cin Bar ai, 


where 
: 5 Bw AP 


f. di y" 0. « _ Betts ¢ — ee f, ?. 
By -21= + a iF B. Be ? Te @2 


Substituting into equations (13a) and (14a) the coefficients of equation 
(22) become: 
A, =) ley. b, 
AG mae B, _ b, Ba 
A,= 20BR,B, +c B. —b,B, 
A = | — 6,8, 
A = aie & c B, +d = be Ba 





FS 


F. Temperature and Pressure Ratios 


The differential equation for the temperature rise through the 


compressor is integrated to determine the temperature ratios as 


follows: aT < ue ¥ (6) 
dé Cp 
ee 
U 
Tig) — To) = oe T(o) dé 


Substituting the expansion of YT (6) and dividing by T(o) = G® the 


temperature ratio is: 


§ § » pS 
ince U" é ty! _¢t. | Sc 
Ta 7 1 t geste eri ftde seal as eB fei. | 





= £ the temperature ratio becomes: 


Let € = 1 and substitute 
ome =) B 


The flow coefficient is defined as: 


,. Yo. m a mR 
v eUA PUA 





The flow coefficient ratio at any station in the compressor becomes: 


>) _  mRTE) | pw Ue) Av) TE) _pe) Aw) 
bd} mR ps UeAsy) To) pe) Ad) 








since m and U are constant through the compressor. The pressure 


ratio is thus given by: 


(£) T(E) $60) A (ec) (2.4) 


pc) T fo) = O°§) A C8) 








sty, 
The temperature ratio is given by equation (23) and the flow coeffi- 


cients by equations (20) and (22). The area ratio is fixed by the 


geometry at the design condition given by equation (6) as: 


ee ee 


4 g 
JA - -f 4 


Integrating gives: 


6 








A(g) ee 

In A (0) = f, lh Be ile, 
Alo) _ (eee) 
A(s) 6. 


The integrals in equation (23) are found by direct integration 
of by» 6° and o, as given by equations (13) and (14). The results 
are again expressed in terms of the initial conditions and coefficients 


involving the design conditions and axial position (€) through the 


compressor. Thus: 


: j 
f¢, a§ a E P (0) 5 Sc, d (0) 

g y / 2 ‘ 2 
f p d& = a, Plo) + C, P (0) b (0) it d. p (0) 


§ 2 / Zz / 
f > dé = 4G, <0) + C,90 po) + ds D (0) 4 C, ,(0) 
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In terms of the coefficients of equation (12) and the € variable 


coefficients (16) the coefficients of these equations are given by: 





a,— (2 \(Er26,- b,) >, “ane ef (§ a b, } j dm & 


Gay) 2 C-2D_ 5D L 
Q.=— {| ——— ee eee = Es Bos 
: tae { 2 (B i Loti 2) V2: (zB 4,tl em) 


thr - Bast) —Klue- 2B) +(ee8) By 





cw Go (2p _ 1? Ah CQ%+t)  _6D 2d 
tm 7+! Zeti (4,+1)* 2%, (5,41) (x, +1)? 72 











eae +(E +B) = ey } 





D 
ane ZF { € : ear ee | 


Substituting these integrals into equation (23) and weet 
the initial conditions as given by equations (2la) the temperature 
ratio also reduces to a quadratic expression in the velocity perturba- 
tion 6, and the initial flow coefficient perturbation , (0). This 


equation is of the form: 





a | 





Té) 2 2 : 
To) a {(6+6) +D, 6 +D,, 6 + DG, 4,(0+D, $0 +D,9,° | (23a) 


The coefficients of equation (23a) are given by: 








ie E 
D- 7 Be - B. 


Te (4B. +2e,B,,) + te (a, B,') + 








42 





| re ( Be+ g) 





Te. (20,B,B, + B, «}— es 


D,= mu (2 a,B,B, +B, ¢, +B,e}+ +: 





i” (a,B'+B, c+ §) 





D= i (B: a,+Be,+d, +B,e,] fi 


The maximum value of the temperature ratio at a given speed 
(5, = constant) is found by setting the derivative of the temperature 


ratio with respect to , 0) equal to zero: 


( ‘ 
er 3} =P e+ D +20, $i} = 0 


The value of the flow coefficient for the maximum temperature ratio 


is thus: 


(0) ee Dy Git - a Di esmee 
Z 7D, 





The second derivative corresponds to the curvature of the temperature 


ratio curve and is: 





22 





a (eH ee De 
i) (0) To) Ge 


and is independent of B,- The temperature ratio curves can thus be 
completely determined by the coefficients of equation (23a). 

A similar parabolic curve can be found for the pressure ratio 
by substituting equations (20), (22), (23a) and (25) into equation (24), 
and dropping all the terms of By and , (0) higher than the second 
order. The same procedure can be used to express the adiabatic 
efficiency n by substituting the equations of pressure and tempera- 


ture ratio in the equation: 
y-l 
pis) x. | 
pe) 
T (§) 
T (0) 





The compressor efficiency (n.) is determined by setting € = l. 
The second order equations for the pressure ratio and efficiency thus 


reduce to the following form; 


fe : 
pS). (520) {146.6 +E, 0+ E,8da+ E, rE, ge} (25) 


am 7. +O, Pr+ 6,8 +628, giI+C, Gt C,, 0 co» (2.6) 


The coefficients of these equations are most easily expressed in terms 


of the coefficients of the previous equations and are: 




















= Dd. = Ax ° aes |-A, Te D, 
Ra” Et Ge , , >. E+ Be 
As; eerie ineD: 1D 
at a ?. >, (€+8.) iz E+ GB. 
Sam; 2A A.-~ Ax fom Arz ale | D ah Deen); 4.2 | 
Ee >* b, g +B. { 2 Sie ob 











E.- (a) iC 7 - u zee 2 + = anes 


cm (GE Pe Gm ECG e.-2, 


K ny = ah _ Ene. ame Dy, 
c ay {E, = 0) EAD. aa +2D,D,—D 


a 


Z 
K ¥-| _ Wes Sal 
erg | Ss egy) ae 


Cy 
! 





Gre Se Kite —E De £5} 4 DY a), 
in K- | x iT ae | Dy t 


vo (rf) 
— (58 


® 


; 7, = B{K—1) 


The coefficient C, can be assumed to be zero since the 


maximum efficiency at design speed is known to occur at the single 


stage design flow rate. Setting the derivative of the efficiency with 


respect to the initial flow coefficient equal to zero: 


d (oO); = 
ove 14 Ca8, == CG. +r eG: ) O 


the maximum efficiency is found to occur at an initial flow coefficient 


perturbation of: 





dij — { Cas. 


ZCu 


Since the maximum adiabatic efficiency at design speed (6, = 0) must 
occur very close to the design flow, C, should be very small. 
The speed function (B,) for maximum efficiency is found by 


substituting the value of the flow coefficient perturbation ($, (o)) for 


maximum efficiency into equation (26) which becomes: 


— wfi+(c- S84) 6, +(ar Se) @ 


Setting the derivative with respect to By equal to zero the maximum 





efficiency is found to occur at a speed function of: 


= CC zo LZ Cu C2 
?: AC, Cos ane Gta 
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Hl. APPLICATION 


The performance of a ten stage compressor has been computed 
to illustrate the use of the method. The single stage characteristics 
for the research compressor as described in reference 1 have been 
used in this analysis, The single stage characteristics for this 
compressor are shown in Figure 1 in terms of the parameters T 
and f. The values of the design point constants for this example are; 

95, 1 a. 25, Fe =- 37.5, f = 2 Ole mle 
and $, = 0,44. The reference rotor speed for this example is 
Uo = 1000 fps. The gas is air at standard temperature and pressure 
with © = 520°R, C= .24 Btu/Ib°R and ¥% =1.4. The value 
of B and ZB then become; Bo =. OOS mand Zo = - 3.288. The 
broken line curves in Figure 1 are the approximate curves represent- 
| ing the quadratic expansions found by dropping all terms higher than 
the second order in equations (9). This approximation is seen to be 
reasonably accurate for flow coefficients from 0.40 to 0. 50. 

The overall compressor performance is found by letting 
N = 10 and € = 1.0 in equations (23a), (25) and (26). The speed 
variable is computed from equation (19) which can be simplified to 


give: 


a- a+ 42 - a+ (4) 42 - elfal—i] 








The resulting temperature ratio curves for the ten stage compressor 
are shown in Figure 2 as a function of the inlet flow coefficient. The 


pressure ratio and efficiency curves are shown in Figure 3. The 
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maximum efficiency is found to occur at B, = 0. 363 which corre- 
sponds to U/U, = 0. 903. 

The operating point of each stage is found by varying & inten 
equal increments and solving equation (22) for the flow coefficient. 
The operating points at design speed for the first, second, fifth and 
last stage of the ten stage compressor are shown in Figure 4, The 
stage flow coefficient (6) is plotted as a function of the compressor 
inlet flow coefficient (ob (o)). The stage efficiency (n), work coeffi- 
cient (WY) and pressure coefficient (uw) are shown as functions of 
the stage flow coefficient (>). These curves indicate the range in 
which the quadratic approximation is reasonably accurate. The quad- 
ratic expansion of 7 and f are shown in Figure 1 to be reasonably 
accurate for stage flow coefficients between 0.40 and 0.50. The tenth 
stage is seen to be outside this range at design speed for compressor 
inlet flow coefficients greater than 0, 455 and less than 0. 400. 

The last stage flow coefficient will have the greatest variation 
from the compressor inlet flow coefficient for off design speed 
conditions. The operating points for the tenth stage at varying speeds 
are shown in Figure 5, 

The operating range of the tenth stage indicates that the method 
should be reasonably accurate for a fairly wide range of compressor 
inlet flow coefficients at design speed. For varying speeds reasonable 
accuracy can be expected at speeds varying up to ten per cent off 


design with a similar flow coefficient deviation. 
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CONCLUDING REMARKS 


It has been shown that the off design characteristics of a multi- 
stage compressor can be predicted by a method using approximate 
analytical integration. The method can be expected to give reasonably 
accurate results at design speed for compressor inlet flow coefficients 
varying as much as four per cent from the design condition. At off 
design speeds the method should be accurate for speeds varying up to 
ten per cent from design with similar flow coefficient variation for each 
speed. A comparison with actual test results has not been made since 
none could be found for a compressor with all stages having identical 
characteristics. 

This method should be extended to include varying stage 
characteristics. This could be done by assuming a given single stage 
BM ractemietic over a section of the compressor and different single 
stage characteristics over other sections. The approximate integra- 
tions could then be carried out over the various sections using the 


appropriate stage constants for that section. 
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